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ABSTRACT 

We use Canada— France— Hawaii Telescope Legacy Survey data for 170 deg 2 , 
recalibrated and transformed to the Sloan Digital Sky Survey ugri photomet- 
ric system, to study the distribution of near-turnoff main-sequence stars in the 
Galactic halo along four lines of sight to heliocentric distances of ~ 35 kpc. We 
find that the halo stellar number density profile becomes steeper at Galactocen- 
tric distances greater than R ga \ ~ 28 kpc, with the power law index changing 
from n inner = — 2.62±0.04 to n out er = — 3.8±0.1. In particular, we test a series of 
single power law models and find them to be strongly disfavored by the data. The 
parameters for the best-fit Einasto profile are n = 2.2 ± 0.2 and R e = 22.2 ± 0.4 
kpc. We measure the oblateness of the halo to be q = c/a = 0.70±0.01 and detect 
no evidence of it changing across the range of probed distances. The Sagittarius 
stream is detected in the I = 173° and b = —62° direction as an overdensity of 
[Fe/H] ~ —1.5 dex stars at _R ga i ~ 32 kpc, providing a new constraint for the 
Sagittarius stream and dark matter halo models. We also detect the Monoceros 
stream as an overdensity of [Fe/H] > — 1.5 dex stars in the / = 232° and b = 26° 
direction at i? ga i < 25 kpc. In the two sightlines where we do not detect signifi- 
cant substructure, the median metallicity is found to be independent of distance 
within systematic uncertainties ([Fe/H] ~ — 1.5 ± 0.1 dex). 

Subject headings: Galaxy: halo-Galaxy: stellar content-Galaxy: structure stars: 
statistics 
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1. Introduction 

Studies of the Galactic stellar halo set constraints on the formation history of the Milky 
Way and galaxy formation processes in general. For example, contemporary simulations of 
galaxy formation predict that stellar halos of Milky- Way-type galaxies are assembled from in- 
side out, with the majority of mass (50% — 80%) coming from several massive (10 8 — 10 10 M Q ) 
satellites that have merged more than 9 Gyr a go, while the remaining ma s s comes from lower 
mass satellites accreted in the past 5—9 Gyr ( iBullock fc Johnston! 120051 ; iDe Lucia fc Helmi 
20081 ). The actual fraction of massive versus less massive mergers will depend on the forma- 
tion history of the galaxy in question. 

A further prediction is that the stellar halos should be more centrally concentrated 
and should have steeper densit y profiles at moderate rad ii (> 20 — 30 k pc), than host dark 
matt er halos (e.g., F igure 9 in Bulloc k fc Johnston! 120051 or Figure 12 in IDe Lucia &: Helmi 
2008 ). According to Bullock fc Johnston (j2005l ). "the difference in profile shapes — and the 
steep rollover in the light matter at moderate to large radii — is a natural consequence of 
embedding the light matter deep within the dark matter satellites: the satellites' orbits 
can decay significantly before any of the more tightly bound material is lost." Hence, they 
anticipate a correlation between the extent of the stellar halo (steepness of the density profile) 
and the extent (or mass) of satellites that built the halo: less extended (more massive) 
satellites will build more concentrated stellar halos. Therefore, the contribution of massive 
mergers to the formation of the Milky Way halo can be constrained by characterizing the 
stellar halo number density profile over large distances and over a wide sky area. 



The large area covered by the Sloan Digital Sk y Survey (SPSS: iYork et al.ll2000l ). with 
accurate photometric measurements (~ 0.02 mag, llvezic et al.l I2004J ) and faint flux limits 
(r < 22), allowed for a novel approach to studies of the stellar distributi on in the Galaxy . 



Using a photometric parallax relation appropriate for main-sequence stars, lJuric et al.l (12008 



hereafter J08) estimated distances for a large number of stars and directly mapped the 
Galactic stellar number density to heliocentric distances of 20 kpc. They found that the 
halo stellar number density distribution within 20 kpc of the Sun can be fit with a two 
parameter, single power-law ellipsoid model 



p(R, Z) oc [R 2 + (Z/q) 



2nn/2 



where R and Z are the cylindrical galactocentric radius and height above the Galactic plane, 
respectively, n = —2.77 ± 0.2 is the power law index, and q = c/a = 0.64 ± 0.1 is the ratio 
of major axes in the Z and R direction, indicating that the halo is oblate (flattened in the 
Z direction). However, additional data suggest that the J08 single power law halo cannot 
be extrapolated beyond 20 kpc. A kinematic analysis by Carollo et al. (2007, 2010) suggests 
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that the halo consists of two components with different spatial density profiles and median 
metallicities, with the "inner" to "outer" halo transition happening at 15-20 kpc, and the 
density profile becoming shallower beyond th at point. On the o ther hand, the distribution 
of RR Lyrae stars from the SEKBO survey (jKeller et al.ll2008l ). and RR Lyrae and main- 
sequence stars from SDSS stripe 82 data seem to show a steeper density profile be yond 30 kpc 



(main sequence stars can be detected up to 40 kpc in coadded SDSS stripe 82; ISesar et al. 
2010h . 



While these studies indicate a change in the halo density profile, each has its shortcom- 
ings: kinematic studies do not give a direct measurement of the density profile but rather 
model it under an assumed dark matter halo potential, RR Lyrae star s are relatively sp arse 
tracers of stellar number density (~ 5 kpc -3 in the solar neighborhood; ISesar et al.ll2010i and 
references therein), and the SDSS stripe 82 region covers only about 1% of the sky. Ideally, 
the halo stellar number density distribution should be mapped using main-sequence stars 
over a large fraction of the sky and to a distance of at least 50 kpc. Such st udies will be en - 
abled by next g eneration wide-are a surveys such as the Dark Energy Survey (ILin et al.1120 09'). 



Pan-STARRS (Kaiser et a 



2008b: LSST Science Book 



2009|). 



2002) and Large Synoptic Survey Telescope (LSST; llvezic et al. 



Meanwhile, the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) has ob- 
served 170 deg 2 of sky in four fields as part of CFHTLS "wide" surve}o The names, po- 
sitions, and sky coverage of CFHTLS "Wide" fields are listed in Table [TJ Due to their 
relatively small sky coverage, these fields are effectively "pencil-beam" surveys for the pur- 
poses of this paper. Despite the small sky coverage (comparable to SDSS stripe 82), these 
data are very useful b ecause of their depth (95% completeness at %' = 23.5 for point sources; 
Goranova et al.ll2009l ). corresponding to a distance limit of ~ 35 kpc for main-sequence stars, 
and because they observe lines of sight unexplored by other surveys. These properties allow 
one to study how the halo density profile changes as a function of distance and line of sight, 
both in the inner and outer halo. 

This paper is organized as follows. In Section [21 we give an overview of CFHTLS data 
and describe the synthesis of SDSS ugri magnitudes from recalibrated CFHTLS u*g'r'i'i 2 
observations. Th e recalibration and tr ansformation of CFHTLS data into the SDSS pho- 
tometric system (IFukugita et al.l Il996l ) is done to allow the usage of CFHTLS data with 
relations already defi ned on the SPSS system, such as the color-luminosity and photometric 
metallicity relations (llvezic et al.ll2008al : iBond et al.l 120101 ) . The CFHTLS z'-band observa- 
tions are not recalibrated as they were not publicly available for all CFHTLS fields at the 



1 http : //terapix . iap . f r/cplt/oldSite/Descart/summarycf htlswide .html 
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time of writing. In Section |3l we analyze the distribution of stellar counts as a function 
of position and compare it to the J08 halo model. The best-fit broken power-law model 
is derived in Section HJ In Section we study the metallicity distribution in the halo and 
analyze the dependence of best-fit model parameters on the adopted metallicity distribution. 
We finish by discussing our conclusions in Section |6j 



2. The Data 

We perform several data quality tests and post-processing steps before using the CFHTLS 
data in scientific analysis. First, we photometrically recalibrate the data by utilizing repeated 
CFHTLS observations and then transform them to the SDSS photometric system. We also 
investigate the performance of different types of CFHTLS magnitudes. 



2.1. Overview of CFHTLS Data 



We use the CFHTLS data processed by t he MeqaPipe image p rocessing pipeline (IGwyn 



20081 ). The pip eline takes as input MegaCam ( iBoulade et al.ll2003l ) images detrended by the 
Elixir pipeline (IMagnier fc Cuillandrdl2004j ). and performs an astrometric and photometric 
calibration on them. The calibrated images are resampled and combined into imag e stacks. 
The catalogs of sources are derived by running SExtractor (IBertin fc Arnoutslll996l ) on each 
image stack. The resulting catalogs only pertain to a single band; no multi-band catalogs 
are generated by the MegaPipe pipeline. The astrometry i n these catalogs is accurate to 
within 0".2 relative to external reference frames (lGwynll2008l ). The catalogs are available for 
each pointing (a ~ 1 deg 2 part of a CFHTLS field) through the Canadian Astronomy Data 
Centre (CADC) Web sitefl 

Since multi-band catalogs are not available, we generate them for each pointing by 
positionally matching sources detected in u*g'r' bands to sources detected in the %' or i' 2 
band using a 1".5 matching radius. The %' and i' 2 bands are the deepest of the CFHT bands, 
so most objects detected in other bands will also be detected in these two. The i' 2 filter is 
the new CFHT filter that was installed after the %' filter broke in 2007 Octobejfl We further 
match ~ 1.7 million C FHTLS sources overlapping the SDSS footprint to SDSS DR7 data 
(lAbazajian et al.ll2009f ) using the 1".5 matching radius. 



• fottp : //wwwl . cad.c-ccda.hia-iha.nrc-cnrc . g c . ca/community/CFHTLS-SG/docs/cf htlswide . html 



See http : //wwwl . cadc-ccda. hia-iha.nrc-cnrc .gc . ca/megapipe/docs/f ilters .html 
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2.2. Star-Galaxy Separation 



We separate point-like and extended sources using the half-light-radius (HLR) measured 
(in pixel units) for each detected source by the SExtractor. For poi nt-like sources, the HLR 
is independent of magnitude and depends only on image seeing (jSchultheis et al.l 120061 ). 
For each CFHTLS band and pointing, we remove the dependence of HLR on image seeing 
by subtracting the median HLR value of bright (17.5—18.5 mag) sources. Bright sources 
ar e used in thi s procedure because they are dominated by point-like sources (see Figure 9 



Gwynl 120081 ). We find that after subtraction, the distribution of HLR values of bright 
sources in the r' and i' (or i' 2 ) bands can be modeled as a ~ 0.1 pixel wide Gaussian centered 
at zero. 

We classify a source as a star if its HLR in the r' and i' (or i' 2 ) bands is less than 
0.2 pixels (after removing the dependence of HLR on image seeing). Sources that do not 
satisfy this condition are classified as galaxies. To esti mate the quality of t his classification, 
we compare it to the SDSS star-galaxy classification ( ILupton et al.l 120021 ) based on deep, 
co-added SDSS stripe 82 data. The star-galaxy separation in co-added SDSS stripe 82 data 
is reliable to at least r ~ 23 (Annis, J. et al.(2011), in preparation), and for purposes of this 
comparison we consider it to be the ground truth. 

Figure [T] shows the fraction of SDSS stars identified as stars in CFHTLS data (complete- 
ness) and the fraction of CFHTLS stars identified as galaxies by the SDSS (contamination) 
as a function of r' magnitude. Using this plot, we estimate that the observed number counts 
will be underestimated by about 5% for r' < 21 and overestimated by about 15%— 20% at 
the faint end (r' ~ 22.5). 



2.3. 



Recalibration and Transformation of CFHTLS Data to the SDSS ugriz 

Photometric System 



Following iPadmanabhan et al.l (120081 ). the photometric recalibration of CFHTLS data is 
separated into "relative" and "absolute" calibration. This process is schematically illustrated 
in Figure |2J For relative calibration, we use overlaps between pointings to recalibrate all 
pointings within a CFHTLS field onto a single photometric system (specific to that field). 
This calibration is then tied to the SDSS photometric system (absolute calibration). The 
final step in the process is the synthesis of SDSS ugri photometry from recalibrated CFHTLS 
u*g'r'i'i' 2 observations. 



-6 - 



2.3.1. Flux Extraction 

The CFHTLS catalogs provide adaptive-aperture (m aMio ) and fixed- aperture magnitudes 
(maper) for sources detected by SExtractor. To compare fluxes extracted by these two meth- 
ods, we select stars from the W3 field and calculate m auto — m aper residuals, where m stands 
for u* , g', r', and i' bands. The residuals are binned in m auto bins and median values are 
calculated for each bin. The dependence of median values on m auto for u*g'r'i' bands is 
shown in Figure O In general, we find that the median difference between adaptive- and 
fixed-aperture magnitudes increases linearly towards fainter magnitudes, reaching ~ 0.04 
mag in the i' band. 

This indicates that there is a problem with one or both extraction methods. As we will 
demonstrate in Section 12.41 the adaptive-aperture magnitudes (m auto ) are responsible for 
the dependence seen in Figure [3], based on the following reasons. First, measuring adaptive- 
aperture magnitudes is a more complex process than measuring at the fixed aperture (for 
details see SExtractor v2.5 manuao), making adaptive-aperture magnitudes less robust. 
Second, a caveat in the SExtractor v2.5 manual warns of potential problems with adaptive- 
aperture magnitudes when the SExtractor R m i n parameter is set too low. This caveat states 
that, "when signal to noise is low, it may appear that an erroneously small aperture is taken 
by the algorithm. That is why we have to bound the smallest accessible aperture to R m in" ■ 
Therefore, if R m in was set too low during flux extraction, inadequate (too small) apertures 
may have been used, making m au to fainter and causing m auto — m aper residuals to be biased 
towards positive values. Even though the caveat states that this may happen at fainter 
magnitudes (low signal-to-noise), Figure [3] seems to indicate that brighter magnitudes may 
be affected as well. We therefore use fixed-aperture magnitudes in the rest of this work. 



2.3.2. Relative Calibration 

Relative calibration places all pointings within a CFHTLS field onto a single photometric 
system (specific to that field) using repeatedly observed stars that are located in regions 
where pointings overlap (CFHTLS pointings overlap about 3" in right ascension and 4" in 
declination directions). Mathematically, the problem of relative calibration can be described 



4 https : //www. astromatic .net/pubsvn/sof tware/sextractor/trunk/doc/ sextractor.pdf 



-7- 



as a x 2 minimization problem ( IPadmanabhan et al.l 120081 ) with 



rista 



X ^ ] Xi ^ ] ^ ] 

i i jeO(i) L 



m j m j,orig ^j,f 



(2) 



where n star is the number of unique, repeatedly observed stars in the / = Wl, W2, W3, 
W4 field, j runs over multiple observations (pointings), 0(i), of the ith star with unknown 



true magnitude m 8 



m 



orig 



and a are the magnitude and its error supplied by MegaPipe, 



respectively, and A^ e j is the correction needed to place the jth pointing on the / field's 
relative photometric system. The system defined with Equation [2] is overdetermined, since 
the number of unknowns (n star + n [parameters] ) is smaller than the number of observations 
(n b s is at least 2n star ). Equation [2] can be e xpressed in matrix form an d can be solved using 
sparse matrix techniques (see Section 3.2 in IPadmanabhan et al.ll2008l ). 



To solve Equation [2] for each field, we use repeatedly observed stars with 17 < m orig < 23, 
where m orig = u*, g', r', i', i' 2 . The observations in this magnitude range are not saturated and 
their reported photometric errors are smaller than 0.2 mag. The observa tions are supplied 
to a sparse matrix inversion code written by IPadmanabhan et al.l (120081 ). and the best-fit 
A r j e j values returned by this code are renormalized so that the median value of A^ e j is equal 
to zero. On average, we find the rms scatter of A^ e l values to be between 0.02 and 0.03 mag, 
reflecting the systematic uncertainty in CFHTLS magnitudes. 



2.3.3. Absolute Calibration 

Following relative calibration, we tie relative photometric systems to the SDSS photo- 
metric system using CFHTLS stars matched to SDSS stars. The goal is to find the best-fit 
A'f s value for each field, such that 

1 2 

m i>re i - m i>sdss - C color - Af s 

/ = (3) 

/ 2 2 
A/ °i + a i,sdss 

is minimized, where i runs over CFHTLS stars in the / field matched to SDSS stars, m re i is 
the CFHTLS magnitude after relative calibration (m re ; = m or i g + A^ e j), and m s d ss and a s d ss 
are SDSS point-spread function (PSF) magnitude and its error, respectively. For m re i = 
u*, g', r', i', i' 2 , the SDSS color and band are color= u — g,g — r,g — r,r — i,r — i and m s d ss = 
u,g,r,i,i, respectively (colors and magnitudes are not corrected for interstellar medium 
(ISM) extinction). 



Wi 

* 2 =££ 

f=Wl i 
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The C color term corrects the linear dependence of rrii jre i — mi )S d ss residuals on color 
caused by d ifferences between CFHTLS and SDSS spectral response curves (see Figure 1 in 



Gwynl 120081 for a comparison of spectral response curves). This term reduces the scatter in 
m i,rei— m i,sdss residuals and improves our estimate of Aj bs . Since CFHTLS and SDSS spectral 
response curves should not change significantly with time, the C color term is simultaneously 
fit for all four CFHTLS fields. Note that this term is only used when estimating A^; linear 
and other higher-order color terms that model the transformation of CFHTLS bands into 
SDSS bands are derived in Section [2.3.41 

To determine A^ values, we use CFHTLS stars with m re i > 17 matched to SDSS stars 
with irtsdss < 21 and with the SDSS color as specified in Table [2j The best-fit C and A^ 5 
are also listed in Table [2J Finally, we define recalibrated CFHTLS magnitudes as 

mcfht = m jJt0ng + A$ + Af s = m jJ>orig + Afj, (4) 

where A*°| is the total correction for the jth pointing in the / field. The A l °j values for 
CFHTLS fields are listed in Tables El 



2.3.4- Transformation to SDSS Bandpass es 

With CFHTLS m c fht = u*, g', r', i', i' 2 observations calibrated onto the SDSS system, we 
now derive the equations that transform CFHTLS observations into SDSS m s d ss = u, g, t, i, i 
magnitudes. In general, the transformation from CFHTLS to SDSS bandpasses can be 
defined as 

m = m cfht + /(color) + Z , (5) 

where Zq is the constant term, and f (color) is some function of CFHT colors (colors and 
magnitudes are not corrected for ISM extinction). 

To find f (color) and Z Q for each SDSS ugri bandpass, we bin m s d ss —rricfht residuals as 
a function of CFHT color, and fit polynomials to m s d ss — m c fht medians. Here we only use 
CFHTLS stars with m c fht > 17 matched to SDSS stars with m s d ss < 20. The dependence 
of m s d ss — m c fh t residuals on CFHT color is shown in Figure 0J 

The best-fit polynomials shown in Figure H] define the transformation of recalibrated 
CFHTLS magnitudes into SDSS magnitudes: 



u 

9 
r 



= u* + 1.07 - 2.375(w* - g') + 1.954(u* - g'f - 0.483(w* - g'f 
= g' + 0.05 - 0.062(0' - r') + 0.365(#' - r'f 
= r' - 0.05 + 0.275(0' - r') - 0.380(#' - r'f 



(6) 
(7) 
(8) 
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i = i' - 0.002 + 0.092(r' - i') - 0.015(r' - i'f (9) 
2 - 0.005 + 0.145(r' - i' 2 ) - 0.280(r - i' 2 f + 0.140(r - i' 2 f . (10) 



i = i 



Using Equations 161 to [T0| we synthesize SDSS photometry (m = u,g,r,i) for all CFHTLS 
observations. Note that these transformations can be used even if the relative recalibration 
step from Section 12.3.21 is skipped because the best-fit A™j values are renormalized so that 
their median is zero. 



2.4. Quality of Photometric Calibration 

To estimate the quality of synthesized photometry, we use CFHT sources brighter than 
m = 17 that are matched to SDSS stars with m s d ss < 20, where m stands for u, g, r, and i. 
For each band m we calculate m — m s d ss residuals, and bin them as a function of color (not 
corrected for ISM extinction, color= u — g,g — r,g — r,r — i for m = u,g, r, i, respectively). 
The absolute value of m — m s d ss medians for m = g,r,i bands is smaller than 0.01 mag 
and shows no dependence on color, for all CFHTLS fields. For Wl, W3, and W4 fields, the 
absolute value of u — u s d ss medians is smaller than 0.02 mag, and shows no dependence on 
u — g color, while for the W2 field the u — u s d ss medians seem to show linear dependence on 
u — g color for u — g < 1.3, as illustrated in Figure \5\ (bottom panel). Since all four fields were 
calibrated using the same procedure, and the transformations from recalibrated CFHTLS 
to SDSS magnitudes were derived using the data from all fields, the u — u s d ss dependence 
on u — g color points to a problem with the original CFHTLS w*-band observations in the 
W2 field. We hypothesize that the u — u s d ss dependence on color may be due to incorrectly 
determined color- dependent air-mass term by Elixir or MegaPipe pipelines. 

The systematic uncertainty of synthesized magnitudes can be determined by comparing 
synthesized photometry of repeatedly observed CFHT sources; such sources can be found in 
regions where pointings overlap. We find that the sys tematic uncer tainty is ~ 0.03 mag (see 



Figure ED , which is consistent with the value cited by iGwynl (120081 ) . 



To quantify the non-linearity of synthesized photometry, we bin m — m s d ss residuals in 
m = u, g,r,i magnitude bins. As shown in Figure [?l the medians of g — g s dss residuals in g 
magnitude bins show no dependence on magnitude, indicating linear behavior of synthesized 
photometry. We have repeated this last test using adaptive-aperture magnitudes and have 
found a magnitude dependence in g — g s d ss residuals similar to the one shown in Figure [3j 
Similar results are also obtained for uri bands. These results point to the dependence shown 
in Figure [3] as being due to the problematic adaptive-aperture flux extraction, and justify 
our choice of fixed-aperture magnitudes. 
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3. Analysis of the Number Density Distribution Profiles 

The data presented in previous sections allow us to measure the stellar number density of 
near main sequence turnoff (MSTO) halo stars, and examine it as a function of position in the 
Galaxy. The sample we obtained from the CFHTLS extends to distances and Galactocentric 
radii (R ga i) nearly a factor of two greater than previous wide- area studies that used main- 
sequence stars (e.g., 5 < R ga i < 15 kpc; J08). T he sample also ove rlaps in range with studies 



based on RR Lyrae stars (5 < R ga i < 110 kpc; ISesar et al.ll2010l ). and it therefore presents 



an opportunity to examine the behavior of the halo density profile in the intermediate range 

(5 < Rgai <30 kpc). 

We begin by selecting a sample of near-MSTO stars with the following criteria: 

0.2<#-r<0.3 (11) 
g > 17 & 17 < r < 22.5 & % > 17 (12) 
5 < D/kpc < 35 (13) 



where magnitudes and colors are corrected for ISM extinction using maps from lSchlegel. Finkbeiner fc Davit 
(119981 ). The g — r color cut serves to select near-MSTO stars, and D is the heliocen t ric dis- 
tance to a star computed using the photometric-parallax relation from llvezic et al.l (l2008al . 
see Equations A2 and A7 in their Appendix A): 

M r = -0.56 + 14.32x - 12.97x 2 + 6.127x 3 - 1.267x 4 + 0.0967x 5 - l.ll[Fe/H] - 0.18[Fe/H] 2 . 

(14) 

The x = g — i color in Equation [14] is computed from the more accurately measured g — r 
color using the stellar locus fit from J08 (see their Figure 8). The [Fe/H is estimated from 
the u — g and g — r color using the photometric metallicity relation from iBond et al.l (120 lOf ) 
(see their Equation Al) 

[Fe/H] = -13.13 + 14.09x + 28.04?/-5.51a;y-5.90x 2 -58.68?/ 2 + 9.14a; 2 ?/-20.61a;?/ 2 + 58.20?/ 3 , 

(15) 

where x = u — g and y = g — r. After these cuts, 13692, 7347, 6505 and 6676 stars are left 
in Wl, W2, W3, and W4 beams, respectively. 

For each pencil beam, we bin the resulting subset in ADM = 0.2 mag wide bins in 
distance modulus^, DM = 51og(D) — 5. We thus obtain the distribution of number counts 
AN (I, b, D), as a function of distance D for each pencil beam. 



5 Binning in bins of equal size in distance modulus (as opposed to distance) results in approximately equal 
number of stars per bin, a consequence of the halo density profile being close to -R~ 3 ; power law. 
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We transform the observed counts to density: 

(1 h n \ AN(k, bj, DM) 
PK L h o h u) .21n(10)£) 3 AnADM 1 ' 

where k, hi and AOj are the field centers and area covered by each beam as listed in Table [TJ 
and p(lt, b iy D) is the number density in stars pc~ 3 . The sightlines sampled by CFHTLS data 
are shown in Figure [U 

In panels of Figure [H we plot the dependence of measured number density on the dis- 
tance from the Galactic center, R ga i for each of the four CFHTLS beams. The measurements 
are marked by symbols with error bars and connected with a solid line for clarity. Overplot- 
ted with a dotted line on each panel is the prediction of the axisymmetric oblate halo model 
of J08. Plotted as open circles, and connected by red line segments, are samples within 
\Z\ < 5 kpc of the Galactic plane. As these may be contaminated by disk stars, we leave 
them out of all further analyses. The residuals of the data for the J08 model are shown in 
Figure [10J 

The density profile observed in Wl beam clearly stands out. While it roughly (within 
~ 20%) agrees with the predictions of J08 at R ga i ~ 15 kpc, beyond that radius the observed 
density begins to exceed the J08 extrapolation, peaking with a factor of ~ 2 excess at 
Rgai ~ 28 kpc, and dropping towards the end of the observed range (R ga i ~ 35 kpc). By 
comparing the location of this overd ensity with the best-fit model of the Sagittarius dwarf 



spheroidal galaxy and its tidal tails (ILaw &: Majewskill2010l ). we conclude that the excess is 



due to the leading an d trailing arm of the Sagittarius stream (jlbata. Gilmore fc Irwinlll994 



Majewski et al.ll2003l ) crossing the Wl beam. While the distribution and metallicity of stars 
in this beam may provide useful new constraints for the study of the stream (see Section [H]), 
this (un)fortunate fact makes the majority of Wl data unusable for the study of the smooth 
halo profile. 

Density profiles in W3 and W4 beams agree within ~ 15% with the predictions of the J08 
halo model in 10 < R ga i /kpc < 25 and 8 < R ga i /kpc < 27 ranges, respectively. The observed 
agreement is nontrivial. First, the directions observed by these beams do not overlap with 
the SDSS data used by J08, and therefore test their model in an entirely different part of the 
halo (especially W4 beam). Second, the CFHTLS data cover a significantly larger distance 
range than the data used by J08, thus validating the extrapolation to ~ 10 kpc greater 
distances (and an order of magnitude change in stellar number density). And finally, while 
J08 did use a small (~ 300 deg 2 ) area in the southern Galactic hemisphere to construct their 
model, their best-fit parameters were largely determined by the > 6000 deg 2 of data from 
the north. The fact that W4 beam (b = —41°. 84) is very well matched by the model puts 
a constraint to any asymmetries between the northern and southern hemispheres to < 15%, 
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out to distances of at least R ga i ~ 25 kpc. This conclusion is also supported to distances of 
D < 18 kpc by the analysis of north versus south SDSS III wide-area imaging (Bonaca, A. 
et al.(2011), in preparation). 

Beyond R ga i ~ 25 kpc, in both W3 and W4 beams, the J08 model overpredicts the 
observed counts. In particular, the J08 model overpredicts the density observed in the final 
two bins of W3 beam by 20% and 40%, respectively. Since we have measured incompleteness 
to be on the order of ~ 5% in almost the entire observed range (see Figured]), and have 
made the bins volume-complete, this turnover cannot be a result of observational bias but an 
indication of a change in halo density profile beyond R ga i > 25 kpc. The overprediction by 
the model may be even greater since the observed CFHTLS stellar counts are overestimated 
by 15-20% at the faint end due to the contamination by galaxies, as shown in Figure [TJ 

The profile exhibited by field W2 is also unusual. As shown in Figure O it is significantly 
steeper than either W3 or W4 profile or the J08 prediction. By itself it is well described by 
a single n = —4.5, q = 0.65 power law. This fit, however, is incompatible with observations 
from the other three directions where a shallower profile closer to n ~ —2.8 is strongly 
preferred. The curious behavior appears to be a combination of two effects: i) the overdensity 
created by the Monoceros stream at R ga i < 25 kpc and ii) the steepening of the halo profile 
beyond R ga i ~ 30 kpc. 



The Monoceros stream (INewberg et al.ll2002l ) is present in the general direction of W2 



field (7 = 2 3 2°, b = 26°). For example, in an AAT/WFI survey of the anticenter region, 



Conn et all (120071 1 detect the stream in the I = 220°, b = 15° direction at D = 11 ± 1.6 
kpc, as well as at Z = 240°, b = 10° with D = 13.8 kpc. Similarly, J08 are able to trace 
the stream in SDSS star count maps to I ~ 230° and importantly, show that it extends to 
at least b ~ 25° of Galactic latitude where it forms a factor of ~ 1.5 — 2 overdensity with 
respect to the extrapolation of the smooth stellar background before and after the stream 
(e.g., as seen in the top left panel of Figure 13 in J08). This is consistent with the factor of 
~ 2 overdensity observed in W2. Secondly, both of these studies estimate the Galactocentric 
distance and width of the stream of R ga i ~ 18 kpc and AR gai ~ 3 — 4 kpc in direction of 
W2. Given that these characteristics are broadly consistent with the enhancement in the W2 
direction for R ga i < 25kpc, and having no additional evidence to the contrary, we interpret 
the observed enhancement as a detection of the Monoceros stream. 
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4. Detection of a Break in the Halo Density Profile 



Taking into account the enhancements due to the Sagittarius and Monoceros streams, 
a single power law remains an appropriate description of the smooth halo component to 
Galactocentric distances of R ga i ~ 25 — 30 kpc. Beyond this limit, however, the observed 
profile appears to turn over rather quickly and the model needs to be modified to explain it. 

To assess the character of the observed turnover, we fit a series of models of increasing 
complexity to the observed data. We use \ 2 P er degree of freedom (xL/) as the goodness 
of fit metric, and search for minima in x 2 hypersurfaces using a Levenberg-Marquardt non- 
linear solver as implemented by the GNU Scientific Librarjo To increase the likelihood of 
finding the true global minimum, we repeat the minimization procedure with 10,000 different 
initial conditions selected randomly from a plausible range of initial values of each parameter. 



We begin by fitting a single J08-type power law to all admissible data pointsj of beams 
W3 and W4, the first eight data points of Wl beam (those that show no contamination by 
the Sagittarius stream), and the last six data points of W2 beam (those that we judge are 
past the influence of the Monoceros stream). We obtain po — 1-7 x 10 -6 pc~ 3 , q = 0.72, 
n = —2.9 as the best, but less satisfactory (xL/ = 6-8), fit. The p is the normalization 
(number of stars per pc 3 ) for the 0.2 < g — r < 0.3 color bin we use. For comparison, the 
fiducial model with J08 parameters shown in Figure M has Xdof = 9.1 when fitted to the 
same data. 

We next increase the complexity of the model by allowing for triaxiality of the ellipsoid, 
parametrized by w = b/a (the ratio of ellipsoid axes): 



This addition makes practically no difference; the best-fit model changes only slightly (po = 
1.7 x 10~ 6 pc -3 , q = 0.72, n = —2.96, w = 1.02) while xlof actually increases (to ~ 7) 
because of the extra degree of freedom. 

We continue by permitting the triaxial halo ellipsoid to rotate in the X — Y (Galactic) 
plane by an angle <fi: 




p{x,y,z) oc (x H + - ) 2 



(17) 



w q 



x' = cos(0) x — sin(0) y 
y' = sin(0) x + cos(0) y 



6 GSL version 1.13; http://www.gnu.org/software/gsl/ 



7 Those having \Z\ > 5kpc to avoid any contamination by the disk. 
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/ 2 2 

p(x,y,z) oc (x' 2 + - + - )* (1* 



This five-parameter model marginally improves the fit (Xdof = 6-8), but converges to param- 
eters po = 1 x 1CT 6 pc~ 3 , q = 1.11, n = —3.3, w = 1.3, = 230° that are strongly excluded 



by prior data (e.g., IChen et al.ll200ll . J08, and others). 



Given the results of this series of experiments, a single power law is unlikely to describe 
the observed counts: the data require a functional form allowing for a change in the radial 
profile beyond R ga i > 25 kpc. We therefore attempt a series of simple "broken power law" 
models, where the density follows one (the "inner") power law until radius i?br is reached, 
and the other (the "outer" power law) beyond. 

We begin with a minimal extension of a single power law model, allowing for a change 
of the power law index beyond a certain radius i?br : 

R e = (x 2 + y 2 + - 2 )^ 
Q 

o(x y z) oc { (i?e)nmnCr ' Re < (19) 
p{x,y,z) oc | (jRe)n _ Re>R ^ (19) 

Note that, as the ellipsoid is allowed to be oblate or prolate, the break radius R\> T is only 
equal to the physical Galactocentric radius R ga i on the x axis (along the line connecting the 
Galactic center and the Sun). In the vertical direction, the physical radius corresponding to 
Rhr is reduced by a factor of q~ x . 

The above model, with five free parameters, produces a significantly better fit to the 
data (Xdof = 3.9). The best-fit parameters for the inner power law, po = (1.45 ±0.05) x 10~ 6 
pc - 3 , q = 0.70 ± 0.01, n inncr = -2.62 ± 0.04, R hr = 27.8 ± 0.8 kpc, n outcr = -3.8 ± 0.1 are 
in excellent agreement with the J08 model, while beyond i?br — 27.8 kpc the best-fit profile 
becomes steeper than the J08 model. 

As shown in Figure [TUJ the observed profiles are better fit by the broken power law than 
by the J08 model, excluding regions with known tidal streams. A fit that entirely excludes 
Wl and W2 fields (not just the regions with known tidal streams) does not strongly constrain 
the broken power law model. In this fit, the normalization p has high (~ 50%) fractional 
uncertainty and strongly correlates with oblateness q and break radius R\, r (correlation co- 
efficient is ~ 1). This strong correlation is caused by similar positions of W3 and W4 beams 
in x — z and y — z planes with respect to the Galactic plane (see Figure E]). In comparison, 
our best-fit profile that uses all four beams while excluding regions with known tidal streams 
shows much weaker correlation between p , q, and R^ T (~ 0.4). 

The inner parts of the W3 field in our best-fit model do show a systematic underestimate 
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of the counts by the model on the order of ~ 15% (at an approximately la-2a leve@), while 
the opposite occurs in inner parts of the Wl field. These lines of sight point toward high 
latitudes in the northern and southern Galactic hemispheres, and the observed difference 
may be a signature of slight north-south asymmetry of halo star counts recently seen in 
SDSS III data (Bonaca, A. et al.(2011), in preparation). 

To further assess the robustness of the detected break, we examined two more variants 
of the model: the first, where we allowed the outer halo to have an oblateness parameter 
q different from that of the inner halo, and the second, where we fixed the parameters of 
the inner halo to J08 values, and allowed only those of the outer halo to vary. Both cases 
resulted in a similar value of Xdof> as we ^ as similar break radii (i?br ~ 28 kpc) and outer 
power law indices (n outcr ~ —3.8). Importantly, the model with varying qm ner and g ou ter 
produced best-fit values of 0.71 and 0.69 for the two, respectively, indicating that there is 
no evidence for a change in oblateness of the halo across the range of distances examined. 

Based on the above analysis, we conclude that the detection of steepening of the halo 
density power law is robust. To explain the CFHTLS data, the power law index needs 
to change from n inner = —2.62 to n outer = —3.8 around i?b r ~ 28 kpc. An R— Z plane 
visualization of the J08 power-law halo model and the broken power-law model presented in 
this paper is shown in Figure [TTJ 



We also fit Einasto's model (lEinastolll965l ) t o our data to allow easy comparison with 



densi t y profiles obtained f r om iV-body simulatio ns (INavarro et al.ll2004j ; iDiemand. Moore &: Stadel 



2004 



Merritt et al.l 120051 ; iGraham et al.ll2006l ). The best-fit parameters for the Einasto's 
model are n = 2.2 ±0.2, R e = 22.2 ±0.4 kpc, p = 1.06 ±0.05 stars pc~ 3 , and q = 0.70 ±0.01 
with chip d f = 4.25. 

Due to contamination by galaxies at the faint end (see Figured]), the n outer = —3.8 ±0.1 
power law index given above is likely somewhat shallower than what it should be. To 
estimate by how much n outer is shallower due to contamination by galaxies, we determine 
f(r) = [1 — contamination(r)]/completeness(r) for each star, where r is a star's r band 
magnitude, and completeness(r) and contamination(r) are the solid and dashed lines from 
Figured! respectively. We sum /(r) values in distance modulus bins and use them instead of 
raw number counts when fitting the model. With this approach, a star located in a distance 
modulus bin with higher galaxy contamination will contribute a value smaller than one 
towards the total count in that bin. Finally, we find that n outeT decreases by 0.1 to n outeT = 
—3.9, which is still within the uncertainty of n outer = —3.8 ±0.1 determined previously. 



3 Note, however, that the error bars of adjacent bins are highly correlated by observational errors. 
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5. Metallicity Distribution in the Halo and its Impact on the Best-fit Model 



In addition to stellar number density distribution, the metallicity distribution can also 
provide strong constraints for halo formation models. Thanks to CFHTLS u band data, 



i t is po ssible to compute photometric metallicity using a method developed by llvezic et al. 



(j2008at ). Furthermore, the absolute magnitude of a star, and hence its distance, depend on 
the star's metallicity (Equation [H]) and thus it is important to study systematic errors in 
best-fit model parameters as a function of metallicity distribution. 



Figure [12] shows the median halo metallicity in four CFHTLS beams as a function of 
distance from the Galactic center, R ga i- In the two beams where significant substructure is 
not detected (W3 and W4 beams), the median metallicity is independent of distance and 
averages to [Fe/H] ~ —1.5 dex with a ran ge of 0.1 dex . This average value is consistent 
with the median halo metallicity measured by llvezic et al.l (l2008a| ) using SDSS F- and G-type 
main sequence stars, and the range is consistent with the systematic uncertainty inherent to 
this photometric metallicity method (~ 0.1 dex, llvezic et al.ll2008al ). 



The average metallicity in the Wl beam is also independent of distance and averages 
[Fe/H] ~ —1.5 dex. This trend is probably a coincidence since the Sagittarius tidal stream 
passes through the beam and the stream's metallicity is not required to be [Fe/H] ~ — 1 .5 



dex everywhere. On the con trary, models and observations (jChou et al.l 120071 ; ISesar et al. 



20101 ; lLaw fc Majewskil |2010| and references therein) suggest that the metallicity exhibits a 
gradient along the Sagittarius tidal stream — we simply happen to observe the stream where 
its metallicity is [Fe/H] ~ —1.5 dex. 

The metallicity in the W2 beam is a bit higher ([Fe/H] ~ —1.3 dex) than in other 
beams, even at large Galactocentric distances (R ga i > 28 kpc) where, according to Figure [TOl 
(top right p a nel), th e contribution of the Monoceros stream stars should be small. Since the 
Ivezic et al.l (j2008al ) photometric metallicity method depends on accurate u band photom- 
etry, anything affecting u band measurements will also affect the photometric metallicity 
estimate. As discussed in Section [231 the u band measurements synthesized from CFHTLS 
data may be impacted by some calibration issues in this beam, so these issues are likely 
responsible for the apparently higher metallicity in the W2 beam (a +0.05 mag systematic 
offset in the u band will increase the metallicity by 0.2 dex). 

Our o bservation that the h alo metallicity is independ ent of Galactocentric distance goes 
contrary to lCarollo et al.l ( 120071 ) and lde Jong et al.l ( 120101 ) (hereafter deJIO) results who con- 
clude that the halo metallicity changes from [Fe/H] ~ —1.6 dex to [Fe/H] ~ —2.2 dex at 
Rgai ~ 15 kpc. To see how our best-fit model varies when the metallicity changes from 
[Fe/H] ~ -1.6 dex for R gal < 15 kpc to [Fe/H] ~ -2.2 dex beyond R gal ~ 15 kpc, we 
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iteratively modify the metallicity distribution in CFHTLS beams to reflect the metallicity 
distributions shown in deJIO Figure 7 (right, hereafter [Fe/ H] deJ i Q (D)) . New metallicities 
are assigned to CFHT stars by interpolating [Fe/H] from [Fe/H] de j 10 (D) using initial helio- 
centric distances, D (calculated from Equations [T41 and IT5|) . Heliocentric distances are then 
recalculated using new metallicity values, and the metallicities are again interpolated from 
[Fe/H]dejw(D))- This process is repeated (about 2-3 times) until the fractional difference 
between distances in consecuti ve steps dips below 0.15 (f ractional distance uncertainty for 



main sequence stars is ~ 15%; ISesar. Ivezic fc Juridl2008l ). The broken power law model is 



fit to modified data once the convergence is achieved. 

The best-fit model parameters obtained for deJIO-like metallicity distribution are listed 
in Table HI In addition, Table H] also lists best-fit model parameters obtained for constant 
metallicity distributions (i.e., all stars are assumed to have fixed metallicity) and for the 
best-fit model presented in Section HI Within a range of plausible metallicity distributions, 
the best-fit model parameters do not seem to vary greatly and average around q ~ 0.7, 
ninner ~ -2.4, n outer ~ -3.8, and R br ~ 25 kpc. 



6. Conclusions and Discussion 

We have recalibrated CFHTLS "wide" survey u*g'r'i'i' 2 observations and transformed 
them to the SDSS ugri photometric system. Using a series of tests, we demonstrated that 
synthesized ugri observations, when compared to SDSS observations, show no dependence 
on color or magnitude. The only exception to this are the u band observations synthesized 
from CFHTLS W2 field u* band observations, which show a slight linear dependence on the 
it — g color. Median photometric error in synthesized ugri photometry is ~ 0.03 mag at the 
bright end and ~ 0.1 mag at ~ 22.5 mag. 

By obtaining the synthesized ugri photometry from CFHTLS observations, we were 
able to use the photometric parallax and metallicity relations which allowed us to study 
the spatial and metallicity distribution of near-turnoff main sequence stars in the Galactic 
halo to heliocentric distances of ~ 35 kpc. We find that the halo number density profile 
becomes steeper at Galactocentric distances greater than R ga i ~ 28 kpc, with the power 
law index changing from rii nner = —2.62 ± 0.04 to n outer = — 3.8 ± 0.1. While the best-fit 
model parameters do change slightly depending on the adopted metallicity distribution (see 
Table HJ), we find that a broken power law model is required for a good fit to the data. 
We also find that the best-fit Rb r value cannot be smaller than Rb r = 22 kpc even for the 
most extreme assumptions about the halo metallicity distribution (i.e., constant metallicity 
at [Fe/H] = -2.0 dex from 5 to 35 kpc). 
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This study p rovides furt 



ously detected bv lSesar et al 



l er evi dence for the steepening of the halo density profile previ- 
( 20101) using main sequence and RR Lyrae stars from the SDSS 



stripe 82, and by iKeller et all ( 120081 ) using RR Lyrae stars from the SEKBO survey. This 
result is consistent with predictions of g alaxy formation simulations which find a steepening 
of the density profi le beyond ~ 30 kpc (IBullock fc Johnston! 120051 ; iDe Lucia fc Helmil 120081 ; 
Zolotov et aJbopgh . 



We see no evidence of change in halo metallicity within the range of probed distances. 
The halo metallicity ranges between [Fe/H] ~ —1.4 dex and [ Fe/H] ~ ~l - 6 dex , and 
averages at [Fe/H] ~ —1.5 dex. This result runs contrary to ICarollo et all (120071 ) and 
de Jong et all (120101 ) studies which report a metallicity of [Fe/H] ~ — 1.6 dex within R ga i < 
15 kpc, and [Fe/H] ~ —2.2 dex beyond. Only the in situ spectroscopic metallicities of 
distant main sequence stars may provide a definitive answer to this discrepancy. With the 
sky density of near-MSTO stars at high Galactic latitudes of about 100 stars deg 2 , the 
multi-object capability over a wide field of view would be well matched to such a program. 

While the total sky coverage of the four CFHTLS beams is slightly smaller than that 
of SDSS stripe 82 (220 deg 2 versus 300 deg 2 ), the CFHTLS beams provide a much stronger 
constraint on the oblateness (Z to R semi-major axis ratio) of the stellar halo because they 
probe very different Galactic lines of sight. We find the oblateness to be q = 0.70 ± 0.01 and 
see no evidence of change across the range of probed distances (5 < R ga i /kpc < 35). This 
res ult is quite consis t ent w ith the oblateness of the dark matter halo, qoM = 0.72, obtained 
by lLaw &: Majewski ( 2010 ) using the positions and kinematics of Sagittarius stream stars. 
However, we find the minor axis of t he ste llar halo to be aligned with the spin vector of the 
Milky Way, while lLaw fc Majewskil ( 120101 ) find the minor axis of the dark matter halo to be 
perpendicular to the spin vector of the Milky Way. 

We have detected the Sagittarius and Monoceros streams as excesses of stars in CFHTLS 
Wl and W2 fields, respect ively. These detec t ions p rovide new constraints on models of these 
streams. For example, the lLaw fc Majewskil ( 20101 ) model of the Sagittarius dwarf spheroidal 
galaxy predicts positions and velocities of Sagittarius stream stars in the CFHTLS Wl beam 
(Figure IT3"j) . The spatial distribution of Sagittarius stars is predicted to be bimodal, with 
a narrow peak at R ga i ~ 20 kpc, and a broader peak at R ga i ~ 35 kpc. The observed 
distribution of Sagittarius stars in this region, overplotted in Figure ITBI (top panel), has only 
one peak at R ga i ~ 35 kpc. Unfortunately, we do not have radial velocity measurements in 
this region, and cannot quantify contributions of particular streams (leading an d trailing) 



to th is peak. Multi-epoch surveys, such as the Palomar Transient Factory (PTF; lLaw et al. 



20091 ) and LSST, will enable robust identification of bright tracers such as RR Lyrae stars, 
which can be used to map the velocity structure in this region. The resulting improvement 
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in models of these streams may then help to further constrain the shape, orientation, and 
the mass of the Galactic dark matter halo. 
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Table 1. Overview of the CFHTLS Wide fields 



CFHTLS Field Name 


a a (deg) 


S a (deg) 


\ b (deg) 


h b (deg) 


Sky Coverage (deg 2 ) 


Wl 


34.5 


-7.0 


173.12 


-61.59 


72 


W2 


134.5 


-3.3 


231.78 


26.04 


25 


W3 


214.4 


54.5 


98.70 


58.47 


49 


W4 




1.3 


63.32 


-41.84 


25 



a Equatorial J2000.0 right ascension and declination. 
b Galactic longitude and latitude. 
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Table 2. Af s Values for CFHTLS "Wide" Fields 



A a f bs 



CFHT Band 


SDSS Color 


Range 


C 


Wl 


W2 


W3 


W4 


u* 


1.1 < u 


-3 


< 1.4 


-0.155 


0.087 


0.099 


0.098 


0.097 


9' 


0.3 < g 


— r 


< 1.2 


-0.161 


-0.017 


-0.044 


-0.037 


-0.043 


r' 


0.4 < g 


— r 


< 1.2 


-0.018 


-0.030 


-0.019 


-0.027 


-0.027 


i' 


0.3 < r 


— i 


< 1.0 


-0.081 


-0.020 


-0.015 


-0.022 


-0.030 


i'2 


0.3 < r 


— i 


< 1.0 


0.020 





n/a 


n/a 
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Table 3. Afj al Values for the CFHTLS W1-W4 fields 



Pointing" 


u* 


<?' 


r' 


i' 


Wl+0-2 


0.106 


-0.034 


-0.034 


-0.048 


Wl+0-3 


0.105 


-0.030 


-0.028 


-0.040 


Wl+0-4* 


0.077 


-0.023 


-0.026 


0.001 



a Pointings with names ending with "*" have i' 2 
instead of i' observations. 

Note. — This table is available in its entirety 
in machine-readable form in the online journal. A 
portion is shown here for guidance regarding its 
form and content. 
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Table 4. Best-fit Model Parameters and their Uncertainties 



Adopted Metallicity 


pg (10- 6 pc- s ) 


q 6 


'thinner 


R br (kpc) 


Tlouter 


Xdof c 


Equation 1151^ 


1.45 ±0.05 


0.70 ±0.01 


-2.62 ±0.04 


27.8 ± 0.8 


-3.8 ±0.1 


3.9 


de Jong et al. f2010) e 


1.25 ±0.04 


0.76 ±0.01 


-2.32 ±0.04 


24.3 ± 0.3 


-6.4 ±0.1 


9.6 


[Fe/H] = -l.0 f 


1.02 ±0.04 


0.72 ±0.01 


-2.32 ±0.04 


28.5 ± 0.6 


-3.8 ±0.1 


5.7 


[Fe/H] = -1.5* 


1.32 ±0.05 


0.73 ±0.01 


-2.47 ±0.05 


24.8 ± 0.7 


-3.86 ±0.08 


4.4 


[Fe/H] = -2.0* 


1.51 ±0.06 


0.72 ±0.01 


-2.50 ±0.06 


22.1 ± 0.6 


-3.81 ±0.06 


4.6 



a Numbcr density of halo stars with 0.2 < g — r < 0.3 at the position of the Sun. 
b Oblateness (Z to R semi-major axis ratio). 
c Goodness of the fit (smaller is better). 

d Our best-fit model. Metallicity calculated for each star using Equation [TJ] 
c Metallicity distribution adopted from lde Jong et al. I J201fj) Figure 7 (right). 
f Metallicity fixed for all stars at a given value. 
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Fig. 1. — The fraction of SDSS stars identified as stars in CFHTLS data (completeness, solid 
line) and the fraction of CFHTLS stars identified as galaxies by the SDSS (contamination, 
dashed line) as a function of CFHTLS r' magnitude (not corrected for ISM extinction). 
Using this plot, we estimate that the observed number counts will be underestimated by 
about 5% for r' < 21 and overestimated by about 15%— 20% at the faint end (r' ~ 22.5). 
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Fig. 2. — Schematic of the calibration process. For relative calibration (top), the pointings 
within a CFHTLS field (9 tiles) are recalibrated to a photometric system specific to that field 
using A r j el offsets (see Section 12.3.21 for details). For absolute calibration (bottom), the field 
is calibrated to the SDSS system (light/green plane) using a Aj bs offset (see Section [2.3.31 for 
details). 
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Fig. 3. — Dependence of m auto — m aper residuals of stars on m auto , where m auto are adaptive- 
aperture magnitudes, m aper are fixed-aperture magnitudes, and m stands for u*, g' , r', and 
i' bands (blue, green, red, and black lines, respectively). 
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Fig. 4. — Dependence of m s d ss — m c fht residuals on CFHT color, where m sdss = u,g,r,i 7 i 
are SDSS PSF magnitudes and m c fht = u *, g', r ', i', i'2 are recalibrated CFHTLS magnitudes. 
The m s dss — m c fht residuals are shown as dots (not all are shown for clarity) and the symbols 
show their median values in color bins. The error bars show errors in medians. The dashed 
lines were obtained by fitting polynomials to m s d ss — n^cfht medians. 



- 30 - 





o.io r 




0.05 r 








0.00- 








-0.05 r 




-0.10 L 



0.8 1.0 1.2 1.4 1.6 1.8 
u-g 





o.io r 


oa 


0.05 : 


c/: 
"O 






o.oo- 








-0.05 : 




-0.10 1 



i i — i — n — r i — i — \ : . \ — v- — i — I — i — i i I — i — ~i — x- 



i . . ■ i 



J ■ ■ I I I I I L I L_ 



0.8 1.0 1.2 1.4 1.6 1.8 
u-g 



Fig. 5. — Dependence of u — u s d ss residuals (dots) on synthetic u — g color (not corrected 
for ISM extinction) for the Wl (top) and W2 (bottom) fields, where u and g are synthetic 
observations derived from recalibrated CFHTLS observations, and u s d ss is the SDSS PSF u 
band magnitude. The symbols show u — u s d ss medians in u — g color bins, and error bars 
show the error in median. To guide the eye, the solid lines show u — u s d ss = ±0.01 mag. For 
the Wl field, the u — u a d ss medians are within 0.02 mag and do not depend on u-g color, 
while for the W2 field the u — u s d ss medians seem to show linear dependence on u — g color 
for u — g < 1.3, indicating a possible problem with CFHTLS u* band observations in the 
W2 field. 
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Fig. 6. — Median photometric error as a function of magnitude for synthetic w-band (solid), 
g-band (dotted), r-band (dashed), and i-band observations (dot-dashed). The median photo- 
metric error was calculated as the rms scatter of iri2 — mi residuals in magnitude bins, where 
mi and ni2 are repeated observations of a star. The systematic uncertainty in synthetic ugri 
magnitudes is ~ 0.03 mag, as indicated by the median photometric error at the bright end 
(magnitudes brighter than ~ 20 mag). 
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Fig. 7. — Dependence of median g — g s d ss residuals on g magnitude, where g are recalibrated 
fixed- (stars) and adaptive-aperture CFHTLS magnitudes (open circles), and g s d ss is the PSF 
magnitude measured by SDSS. The error bars indicate errors in medians. This comparison 
of CFHTLS and SDSS magnitudes shows that the behavior seen in Figure [3] is due to 
incorrectly measured adaptive-aperture magnitudes. Similar results are obtained for u, r, 
and i magnitudes. 
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Fig. 8. — Visualization of the geometry of CFHTLS wide survey beams used in this paper, 
overplotted on isodensity contours of the J08 halo. 
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R g „, (kpc) R gol (kpc) 

Fig. 9. — Stellar number density, measured in four CFHTLS wide-area survey beams as a 
function of distance from the Galactic center, R ga i- Open circles denote the measurements 
within 5 kpc of the Galactic plane, where the contamination by disk stars may be greater 
than 10%. For clarity, the symbols have been connected by solid lines. Overplotted as a 
dashed line is the oblate power law halo model from J08. Its overall normalization has been 
adjusted to fit the W3 and W4 data at R ga i < 25 kpc, as well as Wl data points satisfying 
Rgai < 15 kpc (to avoid contamination by the Sagittarius stream). The vertical line shows 
the J08 distance limit. The excess density at i? ga i > 15 kpc in the Wl field can be associated 
with the Sagittarius stream, while the overdensity at R ga \ < 25 kpc in the W2 beam is 
consistent with the location of the Monoceros stream. 
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Fig. 10. — Data/model — 1 residuals of the J08 model presented in Figure M (dashed line) 
and of the broken power law model discussed in Section H] (solid line). The vertical line 
shows the J08 distance limit. The overdensities in Wl (top left) and W2 (top right) beams 
are due to the Sagittarius and Monoceros streams, respectively. Beyond R ga i ~ 35 kpc, the 
broken power law model provides a much better fit to the data than the single power law 
J08 model (e.g., the J08 model overpredicts the halo stellar number density by > 50% in 
W2 and W3 beams at R ga i > 35 kpc). 
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Fig. 11. — R — Z plane visualization of the J08 power law halo model (left) and the broken 
power law model presented in this paper (right). The color encodes the logarithm of the 
number density of halo stars (stars pc~ 3 ) predicted by the model. Overplotted are the 
densities derived from the analysis of CFHTLS data (beams W2, W3 and W4) presented in 
this paper (Wl beam is not shown because of the strong contamination by the Sagittarius 
stream). Note the marked improvement in data- model agreement for the broken power law 
model. 
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Fig. 12.— Median photometric metallicity (symbols with error bars) measured in four 
CFHTLS wide survey beams as a function of distance from the Galactic center, R ga i- The 
error bars show error in the median and the error bar at (6.5, -1.5) shows the system 



atic u ncertainty in the adopted photometric metallicity method (~ 0.1 dex, llvezic et al. 
2008aj). Within R ga i ~ 30 kpc, the median metallicity is independent of distance and ranges 
from —1.4 < [Fe/H ] < — 1.6. The change in metallicity at R ga i ~ 15 kpc, reported by 



Carollo et all ( 120071 ) and lde Jong et all ()2010l ). is not evident. Apparently higher metallicity 
in the W2 beam ([Fe/H] ~ — 1.3 dex) may be due to u band calibration issues (see the text 
for a discussion). 
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Fig. 13. — Galactocentric distance (top) and Galactocentric rest-frame radial velocity dis- 
tribution (bottom ) of Sa gittarius stream stars in the CFHTLS Wl beam, as predicted by the 



Law fc Majewskil (120101 ) model. The solid and dotted lines denote stars in the first wrap- 



around leading and trailing tidal streams, respectively. The dashed line represents stars in 
the leading arm that has wrapped more than 360° around the Milky Way from Sagittarius 
(second wrap-around). The solid circles show the observed distribution of Sagittarius tidal 
stream stars in the Wl beam obtained from data/model — 1 residuals of the broken power 
law model (multiplied by 10). The bottom panel shows that the contribution of different 
Sagittarius streams to the observed distribution at R ga i ~ 32 kpc, overplotted in the top 
panel, can be more easily quantified in the velocity space, since the majority of Sagittarius 
trailing stream stars are predicted to narrowly distribute around V gsr < —130 km s^ 1 , while 
the leading stream stars have a much broader velocity distribution and V gsr > —100 km s _1 . 



